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Abstract 


In this report we will discuss properties of and models for nuclear transi¬ 
tions in muonic molecules formed via collisions of muonic atoms of hydrogen 
isotopes with light nuclei like D and Be. Their importance for nuclear astro¬ 
physics and nuclear physics at super low energies is emphasized. 


I. INTRODUCTION 

In this report we will mainly consider three particle systems of the type hli pi and hBe pi, 
where h= p, d, or t are isotopes of hydrogen, li = or and Be = ^Be or ^^Be. 
There are strong indications that these particles can form molecular three-body resonances 
which, however, have not yet been observed experimentally. One of these indications is that 
dHe pi resonances were seen in collisions of ppi and dpi atoms with ^He and ^He [jl| . In the li- 
or Be-ca.se there is, of course, a stronger repulsion (proportional to Z, the charge of the heavy 
nucleus) which, however, is partly compensated by an attractive force (proportional to Z^) 
due to the polarizibility of the hpi atoms. Theoretical estimations |^-Q, in fact, show that 
hli jJL and hBe pi resonances do exist. It is however not only the existence, but in particular 
the formation probability of such systems, which is to be investigated. 

Up to now, only mechanisms connected with the transfer of energy to the electronic degrees 
of freedom in the li or Be atom (Auger transitions or resonance excitations of atoms in 
the hnal state 0) have been considered in the literature. However, it is to be expected that 
some other, also resonant, mechanisms exist, leading to the excitation of a (in the Born- 
Oppenheimer sense) ’slow’ degree of freedom. That is, the analog of the Vesman mechanism 
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P] for dt^ formation should play an essential role also in the present case. Indeed, li atoms 
are usually not free in a hydrogen medium, but can form molecular hydrides of li, like UH, 
or £22 molecules. The whole spectrum of excitations of such molecules is displaced into the 
region of a few eV. That means, if in the three-particle systems hli ^ states exist with a few 
eV energy, a resonant formation of the Vesman type is possible. This observation explains, 
why it is so important to perform precise calculations of the spectra of muonic hydrides of 
li with the different isotopes of hydrogen and li. We shall return to this problem later on 
and discuss the possibility of a realization of the cycling process in a mixture of deuterium 
and li. 

If the hli fi or h Be fi resonances are formed, nuclear transitions can occur in them. Gener¬ 
ally, these transitions are expected to be suppressed in comparison with the analogous ones 
in the h^He^ or h^He^ molecules, simply due to their larger size. However, one can hnd 
a few peculiar cases, where the nuclear reactions are of ’long range’ in the nuclear scale of 
distances. There are two reasons for the occurrence of such long-range strong interactions. 
The hrst one is related to the presence of a nuclear resonance in the compound systems, close 
to the threshold energy of some two-body channels and sometimes even in coincidence with 
it. In other words, long-range nuclear transitions in the molecular system can be expected, 
when the heavy constituents of a molecule would have a nuclear resonance at zero relative 
energy. In this case we expect a rather large value for the overlap integral, which charac¬ 
terizes the probability of a transition between the molecular and the nondecreasing nuclear 
resonance wave functions. As an example of such a situation, we consider the molecule 
d'^Ii jjL. In this case, the resonance state of the ^Be{b/2'^) nucleus occurs near the threshold 
energy of the two-body channel d + '^Ii. The second reason for long-range nuclear transitions 
is connected with the occurrence of two closely spaced threshold energies for some two-body 
channels. It is almost evident that, if transitions between such two channels take place, a 
small momentum transfer would be transmitted in the transition process, and this implies 
a long range of the interaction. 

As an example of such a situation we consider the threshold energies of the two-body states 
d -|- ^Be and p + ^Be*{2~^) or p + ^^Be and t + ^Be{g.s.). In the hrst example the difference 
between the threshold energies of the systems d + ^Be and p -\- approximately 

equals 0.04MeH. In the second example it is equal to 0.004MeH. Below, we will estimate 
the probabilities of nuclear transitions in the systems d'^Bep and p^^Bep. 

These systems are of interest also from another points of view. There is a considerable 
lack of information about the nuclear interactions in the eV — keV energy region. Such 
energies, however, are attainable in muonic systems. By lack of information, we mean not 
only that the dynamics of nuclear transitions in this energy interval is not well understood, 
but also that the validity of the usual properties of the strong interactions is not guaranteed. 
Properties of strong interactions as charge symmetry, iso-invariance, the character of P- and 
T-invariance (or its violation) are all established mainly in the MeV region and, up to now, 
have only been extrapolated to the low-energy region. 

A further reason for the interest in such molecules is related to questions of nuclear as¬ 
trophysics. Here we mean nuclear reactions which took place in the process of primordial 
nucleosynthesis just after the Big Bang, and those occuring in those stars, where light ele¬ 
ments are produced. One of the problems here is the abundance of Ti in the Universe. In 
order to discuss this problem reliably, we hrst establish that the abundance of Ti seen in 
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the Universe is indeed due to primordial nucleosynthesis. Evidence for this conjecture was 
found some years ago from the observation of the abundance of in old stars |^. 

In Figure 1 the abundance of is shown as a function of the abundance of Fe. It can be 
compared with the abundance of Mg (see Figure 2), an element which was dehnitely not 
created in primordial nucleosynthesis, but in the stars. 

This comparison indeed indicates that the abundance of nuclei, in contrast to the case of 
Mg, is constant and therefore cannot be attributed to the processes occuring in the stars. The 
modern theory [jT^ of the creation of light elements during the primordial nucleosynthesis, 
predicts the following distributions as functions of the nucleon-to-photon ratio p. 

In this theory one has to use the cross sections for the nuclear reactions in the keV range 
of energies. These energies are not accessible in laboratory experiments, when we want to 
study for instance reactions involving Ti nuclei, such as 


p+'Ii 
and t + a - 


■ 2 a 

+ 'y . 


We can try to get access to these data via the the observation of such transitions in the 
muonic systems p'^Ii p and t '^He p 


p’^Ii p 


^Be /i + 7 
^Be + p 


and 


t^Hep 


p + 'y 
+ p 


where the relative energy of the heavy particles in the molecule is arround a few eV only. 
Usually, astrophysicists extrapolate the astrophysical S-factoiQ, dehned as 


-2TTri(E) 


( 1 ) 


into the region of low energies. How dangerous such a process is, can be seen from the 
example of the extrapolation procedure for the S-factor of the reaction 


d T d 


Me + 7 . 


( 2 ) 


In Figure 4 we show the extrapolated S-factor and the S-factor, fitted to the new data at low 
energies |^. It is seen that the difference in the ’old’ and ’new’ S(0)-values is two orders of 
magnitude. 

Another reaction that we discuss in the context of muonic systems is the radiative capture 
of protons by the ^Be nucleus 


p 


^Be 


®S + 7 


( 3 ) 


^The S-factor is assumed to be a slow function of the energy, if no resonances in the region 
concerned are present. 
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It is well known [^] that this transition is the last step in the pp-chain in the sun and is the 
bottleneck, which is crucial for the determination of the flux of high energy neutrinos from 
the sun In laboratory measurements of the cross section of this process only the region 


of energies down to IQQkeV is accessible. An extrapolation to the energies at the centre of 
the sun {l.bkeV) has to be made. Here again one of the possibilities to obtain the value 
of the S-factor in the few eV energy region is to search for this transition in the muonic 
molecule p ^Be p. In this case we do not have a threshold resonance in the p + ^Be system, 
but a very peculiar hnal state. The wave function of the ground state of the nucleus has a 
very long tail due to the (on a nuclear scale) extremely small separation energy of the proton 
in which equals to 0.14MeH. Due to this long tail, we again have a process, which occurs 
at large distances and can be favourable for its observation in the muonic system p ^Be p. 

II. THEORETICAL DESCRIPTION OF CHARGE NONSYMMETRIC MUONIC 

MOLECULES 


In this chapter we describe the treatment of the three-body Coulomb problem hZp. Since 
we have two heavy particles and one light particle, it seems that the Born-Oppenheimer ap¬ 
proach will be reasonable to describe this system. Indeed, some estimations of the prop¬ 
erties of the hHe p systems performed by this method (in the one-channel approximation) 
provided a qualitatively correct description of some properties. Of course, the quantitative 
description requires improvements since, within this approximation, it is impossible to re¬ 
produce even the spectra of the usual muonic molecules of hydrogen. Moreover, since we are 
interested in the probabilities of some partial transitions, say the probability of formation, it 
will be more reliable to use a method taking angular momentum conservation into account. 
Another problem within the framework of the Born-Oppenheimer approach is that it is very 
difficult to treat the contribution from the continuum spectra of the two-center problem 
correctly. Thus, to avoid these difficulties we will use another adiabatic formulation, namely 
the hyperspherical adiabatic approach This method is based on an expansion of the 


three-body wave function into the so-called surface functions [T^ . 

Let us start the description of our formalism regarding Figure 5, where the Jacobi coordinates 
are shown. 

In these variables the Hamilton operator of the three-body Coulomb system takes the form 


H = - 

2m 


-Ar 

2M ^ 




+ 


r \R-(3f\ \R + 'yf\' 


(4) 


where 

= M~^ = + {m. + mz)~^, P= ———, 7 = 1-/? (5) 

^ mf, + mz 

and m^, mh and mz denote the masses of the muon, the hydrogen isotop and the heavier 
nucleus (with charge Z), respectively. Introducing dimensionless variables according to 


R 


X = 


y = —, 

a aa 


( 6 ) 
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with a = and a being the Bohr radius of the {Z/i) subsystem, a = (Zme^) the 
Hamiltonian reads 


H = 


2ma^ \ 


/ _ - A- - - - 


+ 


X Z\ay — (3x\ Z\ay +'-)x\ 


(7) 


This suggests to take {2ma?) ^ as our energy unit. Let us dehne the hyperradius p and the 
hyperangle uj via 


x = pcosuj, y = psinuj , (8) 

with 0 < p < oo and 0 < a; < 7r/2. The Hamiltonian (|^) then takes the well-known form 


op^ p^ p 


(9) 


where the grand angular momentum operator A (hi) is given by 


= 


1 
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1 

cos UJ sin UJ H-^-h ^ 


cos UJ sin UJ duj'^ 


cos^ UJ sin^ UJ 4 


( 10 ) 


lx and ly are the orbital angular momentum operators corresponding to the variables x = x/x 
and y = y/y, respectively. For the potential we have 


V(ff) 


-2 


1 1 
coso; Zlof^sinci; —/3a; coso; 


Z\aysm.uj -|- yxcosci;| 


( 11 ) 


The hve angles symbolized by hi = {uj,x,y) together with the hyperradius p provide a 
complete set of variables to describe the positions of all three particles. The Schrodinger 
equation therefore is 


_p-5/2^ 5/2 


p ] 


T(p,fi) = Evl/(p,f2) 


( 12 ) 


Let us now introduce the surface functions mentioned above. We want to use these functions 
as a basis for the solutions T(p, hi) of (|T^. For this purpose we consider the operator 
containing only the angular part of the Hamiltonian, 


p2 p ] 


Bn{p, fi) 


Un{p) H„(p, n) . 


(13) 


Its eigenfunctions i?„(p, fl) are the so-called surface functions, its eigenvalues Un{p) are called 
eigenpotentials. In equation (|l^ the hyperradius p enters only as a parameter. Since the 
surface functions form a complete set on the sphere of constant values of the hyperradius, 
the three-body wave function can be expanded according to 

. (14) 
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The adiabatic approximation consists in the truncation of this sum into a finite number of 
terms. In contrast to the Born-Oppenheimer approach, not the distance between the two 
heavy particles, but the hyperradius (a measure of the average size of the whole system) 
is the ’slow’ variable. Both variables coincide if the mass of the light particle (here m^) 
tends to zero. The two major advantages of our approach, in comparison with the Born- 
Oppenheimer one, are that here the total angular momentum is a good quantum number 
and that the eigenpotentials Un{p) tend to the exact two-body binding energies in the limit 
of p —> oo. Both properties can be derived from equation ([T3|) . 

Instead of solving the coupled system of differential equations, which is obtained by inserting 
(p^) into the Schrodinger equation (|^), we here only treat the so-called extreme adiabatic 
approximation, which is given by the neglect of the derivatives of the surface functions with 
respect to the hyperradius. The remaining radial equation then is 


f'M = Pnip) - E]Up) , 


( 16 ) 


looking exactly as a radial two-body Schrodinger equation, where Un{p) acts as an effective 
potential. Hence the name ’eigenpotential’. 

The main task now consists in the determination of these eigenpotentials. As it can be seen 
from equation (|^), for small values of the hyperradius, A^/p^ dominates over the Coulomb 
potential V/p. Therefore, the surface functions in this region should be proportional to the 
eigenfunctions of A^, which are the well-known hyperspherical harmonics Y[£](r2) 


Y[£](H) = N[c] cos'^w sin^^o; Pfc'“^"’'"^ = \cos2w) ytX{x,y) . (16) 

These functions are characterized by the set of quantum numbers [P] = {k,lx,ly, L, M}. 
They contain the Jacobi polynomial 


p{OL,P) 


(X) 


{x — 1) “(x -1-1) d 

2^k\ dx^ 


l)fc+a(^-l- l)fc+p 


(17) 


and the bispherical harmonics (eigenfunctions of the squared total angular momentum op¬ 
erator L^) 


ytX(i,y)^ E ( 18 ) 

rrixirLy 

with the Clebsch-Gordan coefficients {IxfOxlymy \ LM). The normalization constant Nyc] is 
given by 


Nyc] = 


\ 


2k\ {lx + ly + 2k + 2) T{lx + ly + k + 2) 

T{lx k T{ly -|- /c -|- |) 


(19) 


Large values of the hyperradius and negative energies (the only energies we are interested 
in) correspond to the physical situation where the muon is bound by one of the positive 
charges. Hence, the surface functions in this region look like channel functions 


$[^.](p, H) = p^/2 77„z^(pcosa;) sin'^ y[f (x, y) . 


( 20 ) 
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Here, 7Zni^{x) denotes a hydrogen-like wave function and [m] the set of quantum numbers 
{n, lx, ly, L, M}. The index ’i’ specifies by which of the two nuclei the muon is bound. 

To represent the surface functions in the whole space of p we use the following ansatz 

=Y.Y. «rxK](p) ^ K[C\{p) , (21) 

*=1 [mi] [£] 


which (inserted in equation (p!3D) yields a generalized eigenvalue problem for the determina¬ 
tion of the coefficients an[mi\{.p) and &n[£](p) and the eigenpotentials Un{p)- 
Within this framework we treated the systems dtp, h^Ii p, h'^Ii p and h'^Be p in the states 
with the total angular momentum L = 0. In all calculations we use 120 hyperspherical 
functions (|^). In the case of {dtp) we take 6 channel functions (^) and for the systems 
of higher charge we take 10 channel functions into account. The latter is useful due to the 
higher polarizibility of the {hp) atom by the nucleus of higher charge. Tables I to IV show the 
binding energies of the various muonic molecules using the extreme adiabatic approximation 
(p!5|) in comparison with the results of References 1,0 


III. ESTIMATIONS OF NUCLEAR TRANSITION PROBABILITIES IN 
MUONIC MOLECULES OF BE ISOTOPES 

We start by considering simpler systems with respect to their nuclear structure, e.g. the 
system d '^Be p, to demonstrate the approximations we make and to discuss their reliability 
. First, consider the structure of the spectrum of the ^Be nucleus shown in Figure 6. 

As one can immediately see that is the case of two thresholds close to each other, which 
correspond to the states d + ^Be and p + ®Bfe(2+). The difference between these threshold 
energies is equal to 0.044MeU, a very small value on the nuclear scale. In such a way, if the 
nuclear transition d + ^Be —>■ p + ^Be{2^) takes place, only a small momentum transfer is 
possible and large distances between d and '^Be are involved. 

Now, let us formulate the description of such a transition in the muonic system 

d'^Bep—^ p^Be{2'^) p . (22) 

We will treat both states (initial and hnal) as four-body systems n + p + ^Be -|- p. The 
internal structure of the ^Be nucleus is not important in the range of energy and momentum 
transfer involved. For that reason we will treat the excited state of the ^Be{2~^) nucleus in a 
simple two-body n + '^Be model. Let us introduce the two sets of Jacobi coordinates, shown 
in Figure 7, which are appropriate for the description of the initial and hnal states. 

The total Hamiltonian can be written in the form 

H = Ho + V^ + V^ , (23) 

where the sum of the kinetic energies Hq is given by 

Hq = ho(r) -|- ho(p) + hQ{R) , (24) 

and the potentials of the Coulomb, U^, and strong, , interactions are 
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v° = + and ^ Vl + V^,^ +, (25) 

respectively. To continue, it is useful to introduce the auxiliary channel Hamiltonians Hi 2 , 
dehned via 


Hi — ho(ri) + ho(pi) + hn^Ri) + V^pifi) 


(26) 


and 


H 2 — ho{f 2 ) + ho(p2) + ho{R2) + 

+ ^»Beti + + ^p^Be 5 

and the corresponding Schrodinger equations for the channel eigenfunction^, 


Hitl^i = EillJi , 


f = 1,2 . 


(27) 


(28) 


As one can see from the equations {W) and (^) the channel Hamiltonians do not contain 
the strong interaction between the clusters d and in the initial state and between p and 
^Be in the hnal one. Therefore they describe the three-body Coulomb systems d ^Be p and 
p ^Be p and the internal motions in the deuteron and ^Be nucleus, and therefore are still 
four-body Hamiltonians. 

To solve the Schrodinger equation 


HT = ET 


(29) 


with the full Hamiltonian (|23|) we use the ansatz 


T = CiiJi + C2i’2 ■ 


(30) 


This ansatz reminds us of the well-known LCAO (Linear Combination of Atomic Orbitals) 
approximation, where the role of the ’atomic orbitals’ is played by the nuclear states of d and 
®Hfe(2+). The approximation ( PPP therefore implies that the internal motion in the deuteron 
does not disturb the internal motion in the ^Be nucleus much and vice versa. Let us discuss 
the reliability of that ansatz for the solution of equation (^). We should hrst note that in 
the usual molecular ion of the hydrogen where the LCAO approximation is qualitatively 
a good one, the size of the atomic orbital has the same order of magnitude as the size of the 
H 2 ion itself, namely tq ~ 10“®cm. In the case of the muonic system the ’orbitals’ {d and 
^Be{2~^)) are of a nuclear size (~ 10“^^cm) and the d’^Bep molecule has a size of ~ 10~^cm 
|T5| . So the orbitals in our case are four orders of magnitude smaller than the molecule and 


we should expect much more reliable results from approximation (|30D than in the case of 
the H 2 ion. To continue the analogy with the H 2 ion we consider the picture, in which the 


^For reasons of simplicity we will neglect for the moment the contributions from the non-point-like 
charge distributions in the d and ^Be nuclei. 
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neutron can be bound in some potential well by the proton (position a) or almost with the 
same binding energy by the ^Be nucleus (position b), as is shown in Figure 8. 

It is well-known in quantum mechanics that in such systems, where the possibility exists 
that a particle can be present in either of the two potential wells, splitting of levels should 
occur. This does indeed happen to a signihcant extent in this case as we shall now proceed 
to demonstrate. From the structure of the channel Hamiltonians Hi ^2 we immediately hnd 
for the eigenfunctions " 01^2 that 

V^i = <F,(ri) ^T\pi.Ri) (31) 

and 

^2 = <h8a(r2) <Fr'(p2,4) . (32) 

Here and $ 85 ^ are wave functions, which describe the internal motion in the d and 
nuclei, while the molecular functions for d'^Be^ and p^Befi are given by and 
^moz, respectively. 

Now, after diagonalization of the total four-body Hamiltonian H with the basis functions 
(0) and (|3^) , we easily hnd the coefficients Ci and C 2 in equation (^0]) and the transition 
matrix element M, 


M ~ (I '02 ) • (33) 

For the reaction rate P we then have 

F = «|(v1/|02)P, (34) 


where k is the muonic molecule frequency. It is easy to see that we obtain two solutions of 
equation (|25|), corresponding to the split levels 


= Ar±(A±0i-f02) , 

where is the normalization constant and 

F±/-Fri2 

^ - Hn-E^ ’ 


(36) 


(36) 


with the matrix elements 


/ = ( 01 I -02 ) 


and 


Hij ^ , 


(37) 


which are indeed multidimensional integrals. For the transition rate ([3^) we obtain 

11 + rA^P 


= K 


l + 2Fe(FA±) + |A±|2 


(38) 


From the structure of the channel functions (pT]) and 
integral I we have the relation 


I) it is easy to see that for the overlap 
J = (0i|02)~y'dp |4>”^°'(p,F = 0)|2 =^3, (39) 
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where the value of S 3 can be interpreted as a probability for the heavy particles in muonic 
molecules to be at small distances (zero distance in this case). It is easy to see that the 
matrix elements of the total four-body Hamiltonian (p3D also can be expressed through the 
overlap integral For example for Hu we have 


H 


11 


+ >S'3(V) 


11 

n Be 


+ 


(40) 


where the are the matrix elements of the strong NBe potentials over the wave function 
of the deuteron <1)^. 

Since the nuclear Vnbe potentials have been determined by experimental! data [^, the most 
uncertain value, dehning the rate of the nuclear transition, is S 3 , the probability for 
the heavy particles to be close to each other. To make a rough estimation we will proceed 
in the following way. On a qualitative level the three-body muonic states like pBe can be 
considered as two-body (p -|- Be) systems, moving in an effective potential produced by the 
muon. In that potential we can calculate the probability for the heavy particles to be at 
small distances. Let us call this value 5*2. Then we can suppose that 


S3^S2 = mR = 0)W 


(41) 


The value of the wave function at small distances can be estimated 


23 


as 


I<l)(fl = 0)1^ < (m . 


(42) 


Here m is the reduced mass of the two particles and AE the separation between the s 
and p levels in the effective potential, which produces the wave function ^{R)- Using this 
estimation for the rate of the nuclear transition P~, the only one available in the process 
d fi —p ^Be p, one obtains M 


P = 4.6 * lO^^sec 


-1 


(43) 


The inequality (l4^) probably overestimates the real value, since the repulsion in the effective 
potential is not fully specihed in Reference p^. 


Now let us come to the description of the nuclear transition in the molecule p ^^Be p ||24 

p^^Bep — >t^Bep + Q, Q ^ AkeV . (44) 


The incredibly small energy release in this transition is related to the peculiar properties 
of the ^^Be nucleus. Indeed, as one can see from Figure 9, a double coincidence takes place 
here: the coincidence of the position of the excited state 8.4774MeU with the threshold 
energy in the channel ^Be + 2 n, and a coincidence of the absolute value of the energy of this 
excited state with the binding energy of the triton (8.48MeU). 

An additional peculiarity of that transition, in comparison with the previous one, can be 
seen from the spectrum of the ^^B nucleus. It turns out that this nucleus has excited states 
with energies which coincide with the threshold energy in the channel p+ ^^Be. Keeping this 
in mind, two essentially different mechanisms of the transition can be formulated, as shown 
in Figures 10 and 11. 
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However, since the spins of the corresponding excited states of the nucleus are rather 
high (9/2’'',5/2“), we will neglect contributions from that mechanism. In other words, our 
transition amplitude does not contain the resonant term from Figure 11, due to a sufficiently 
small angular momentum in the initial state. So, we can proceed in complete analogy with 
the previous case. The only complications which arise, are due to the more complex wave 
function of the nuclear state. To describe the process in the frame of the ’LCAO’ ansatz, 
we have used a four particle {2a + 2 n) model for the description of the ground state of the 
nucleus, and a 2 q! model for ^Be{g.s.). 

All we need now are Na and NN potentials, wave functions of the nuclei ^^Be{g.s.), ^Be{g.s.) 
and the triton and molecular wave functions at small distances. The nuclear wave functions 
were written in a simple form, reproducing only the periphery of the nuclei j^J^. For 
the Na potentials three different but phaseshift equivalent potentials p7|-p9| were applied. 
The results are sensitive to the differences in their shape. The molecular wave function was 
constructed by means of the potential 


V{R) = 4.5 H [— arctan(ai?) — 1] e 

71 


-bR 


IH, 


(45) 


with parameters fitted to the eigenpotential calculated by means of the method of surface 
functions |^. Here, H = me'^/T l = 5626.5eH is the Hartree energy unit for the muon. The 
potential (^), of course, produces a more reliable wave function of the molecule at small 
distances as compared to the estimation of equation As NN potentials the Malfliet- 


Tjon potentials have been used. The results for the reaction rates are shown in Table 
V. The fact that the lowest reaction rate occurs with the Na potential given in Reference 
^ , most likely is a reflection of the repulsive character of that potential. 


IV. NEW CYCLING PROCESSES AND OTHER SPECULATIONS 

Considering the properties of molecular systems with isotopes of li nuclei, we can hope that, 
in principle, new cycles analogous to the one occuring in the dt mixture, can take place. To 
see this analogy more clearly, let us remind ourselves of the picture of the cycle with the dtg, 
molecule, as shown in Figure 12. 

Roughly speaking, the two following parameters of that cycle are specihc only for this system: 
1) the probability of the resonant formation of the dtg, excited stateQ, 2) the probability of a 
nuclear transition in the dtg, molecule which also has a resonance behavior. As it was pointed 
out in the introduction, the possibility of resonance formation of d li systems cannot be 
excluded in collisions of the type d/i + liD. If excited states in dli g, systems within the 
energy range E < 2.5eU exist, the Vesman mechanism can be invoked. {E k, 2.5eU is the 
dissociation energy of the molecule liD from the A'TN state |32|). 

Considering the nuclear resonances we see for both isotopes of li a situation, which might 
be even more favourable for nuclear transitions than in the dig case. Indeed, let us look on 
the spectra of the nuclei ^Be and ^Be, schematically shown in Figures 13 and 14, respectively. 


^Resonance molecular formation, as it is well-known |31|, also takes place in dg + D 2 collisions. 
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We see that in both cases there are threshold resonant states in the compound nuclei 
and '^Be and the positions of these resonances are closer to the corresponding thresholds than 
the position of the ^He{?)/2'^) resonance to the d + t threshold. In the last case the ’center’ 
(on the energy scale) of the resonance state ^He{3/2~^) is 50keV higher than the threshold 
energy. A nuclear transition, therefore, only takes place through the tail of the resonance 
curve. Keeping in mind the above pecularities of the molecular and nuclear structures of 
the systems dfi + we can imagine the cycle in the d + li mixture, shown in Figure 15. 
The third parameter crucial for the occurence of the cycling process is the sticking coefficient 
for the muon in the hnal state. Unfortunately, none of these three parameters have till now 
reliably estimated. 


Other interesting phenomena can be observed in the above systems. We are talking about 
the so-called Zeldovich phenomenon or rearrangement of the spectra of Coulomb systems 
by a short range interaction. As it was shown by Zeldovich |Q, in systems with long 


and short range interactions a rearrangement of the spectrum of the purely long range 
potential takes place, if the short range interaction has a resonance or bound state near zero 
energy. In this case the strong interaction cannot be treated perturbatively despite its short 
range character. Both muonic systems d^Ii and d’^Ii ^ are indeed good candidates for 
the manifestation of the Zeldovich phenomenon, since the short range nuclear interactions 
must generate a resonance behavior at zero relative energy of d and li. It is necessary to 
emphasize that the rearrangement of the spectra takes place due to the large size of the 
hadronic subsystem at the conditions we have described above. 

Considering the nonperturbative situation concerning the influence of the strong interac¬ 
tion, we can look for the possibility of nuclear transitions in the corresponding electronic 
molecules. As an example let us consider a molecule of ordinary light water, H 2 An 
interesting property of that molecule is the coincidence of its energy with the energy of the 
excited state of the nucleus ^We. Let us look at the spectrum of this nucleus, schematically 
shown in Figure 16. 

As one can see, the nucleus indeed has a state with an energy incredibly close (within 
the accuracy of all measurable hgures |^) to the treshold energy for the three-body channel 
^^+2p. Since the binding energy of water is only a few eU, we claim that the water molecule 
in the rotational state 1“ and the nucleus ^®Ae in the excited state (1“) of 4.522MeU are 
degenerate states of the same Hamiltonian, describing 18 nucleons and 10 electrons. This 
implies that the wave function of molecular water is not a pure state, but always contains 
an admixture of the nuclear state. Since this excited state of the nucleus can decay to 
other channels, for example, in the channel -|-p, molecular water should slowly disappear 
through the resonance state of the nucleus ^We. In other words, the smouldering (slow 
burning) of water should take place with an intensity dehned by the overlap of the molecular 
wave function and the nondecreasing wave function of the resonance nuclear state. Exactly 
the same is valid for the molecules of hydrides of Li due to the coincidence of the nuclear 
resonance and threshold energies. 
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V. CONCLUSIONS AND OUTLOOK 


In this paper we reviewed work performed by our group during the last few years, devoted 
to charge nonsymmetric muonic molecules. For their treatment we have to consider both 
Coulombic and nuclear few-body problems. We have given special attention to examples 
where the interplay between the Coulomb forces and the nuclear forces is of atmost impor¬ 
tance. We have performed calculations for some of these examples which have been reviewed 
here and others will be published elsewhere. 

For the treatment of the purely Coulombic systems, we employed the method of hyperspher- 
ical surface functions. We have calculated eigenpotentials and eigenfunctions for all isotopes 
of hydrogen and Li and some isotopes of Be. Using the extreme adiabatic approximation, 
our calculations yield no loosely bound states (~ 1 or 2eU) having total angular momen¬ 
tum L = 0. The nuclear transitions were calculated in a scheme motivated by the LCAO 
method, but improved by taking into account all degrees of freedom of light and heavy 
particles. To reduce the numerical effort to manageable proportions we introduced effective 
potentials between the nuclear clusters. Obviously, there is a degree of nonuniqueness in 
this procedure. However the main features entering into our calculations, like the positions 
of the thresholds, are not effected by this approximation. 

We have considered nuclear transitions in the molecules d ^Be /i and p ^^Be /x, where the energy 
releases as well as the momentum transfers are extremely small on the nuclear scale. For 
these reasons the nuclear interactions are expected to have a long range character, which is 
important in muonic systems. 

We now briefly discuss the reliability of the nuclear reaction rates performed in these two 
cases. The crucial ingredient in these calculations is the value of the molecular wave func¬ 
tion at the origin. In the first instance we calculated this quantity by using a rather crude 
approximation, where the Coulomb interactions were only taken into account in an indi¬ 
rect way suggested by Ref. ||2^. In the second example the Coulomb repulsion at short 


distances was taken into account rather reliably, since the eigenpotentials obtained via the 
hyperspherical surface function method have the correct behavior near the origin. 

To proceed further we have to keep in mind the difficulties associated with the Coulombic 
three-body problem in which we are interested. In the first place we should emphasize 
that the systems of the type hZp with Z > 3 are highly excited resonances and up to 
now no complete solution of neither the Schrodinger equation nor the Faddeev equations 
has been presented anywhere. The difficulties of the solution of the Schrodinger equation 
are related mainly to the very complicated boundary conditions in the configuration space. 
Faddeev equations for Coulomb forces are not of Fredholm type in this range of energies and 
should therefore be modified to obtain unique solutions. Additional difficulties arise from 
the fact that the calculations of the eigenvalues must be performed to a very high accuracy. 
As indicated above, the dissociation energies of the molecule liD lie around 2.5eU. That 
means that the accuracy of the calculated binding energies of the systems dli p should be 
of the order of leU or even better. 

It seems that at the present time only two approaches are sufficiently promising for the 
treatment of such systems. The first one is the hyperspherical surface function expansion, 
briefly presented here. From physical reasons we consider resonance states in the extreme 
adiabatic approximation as real bound states having no width. Indeed, the inclusion of 
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the coupling between the channels corresponds to taking very high excitation energies into 
account. The width appears only, when nonadiabaticity is taken into account and a coupled 
system of equations is employed. 

The second promising method, which is very often used in molecular physics |^, also is 
based on the treatment of the resonance state as a bound state. But the realization of that 
idea is more motivated by mathematical arguments. To get a square integrable solution in 
this approach, we use coordinate rotation in the complex plane. It can be shown [M that 
the rotated Hamiltonian has the same eigenvalues as the physical one. 

Concerning the study of the nuclear interaction in the eV — keV energy region, one can 
say that it is the only region of energies in the held of nuclear physics, where there are 
practically no data available. We emphasize that the research proposed here can shed 
light on the properties and behavior of nuclear systems and interactions in the very far 
periphery at distances, exceeding the usual range of the nuclear interaction by several orders 
of magnitude, if one considers muonic systems. 
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TABLES 


TABLE I. Binding energies of {dtfi) with L = 0 

state 

ground state 

excited state 

present 

-295.54 

-35.80 

H 

-319.14 

-34.83 

TABLE 11. Binding energies of {h^Ii fi) with L = 0 

system 

p^Ii H d^Ii p 

t^Ii p 

present 

-24.3 -23.8 

-35.3 

i 

-17.6 -18.5 

-19.8 

TABLE III. Binding energies of {h fi) with L = 0 

system 

p’^Ii H d'^Ii p 

p 

present 

-20.8 -25.9 

-37.5 

i 

-21.0 -22.0 

-23.3 

TABLE IV. Binding energies of [h ^Be fi) with L = 0 

system 

p^Be p 

d'^Be p 

present 

-11.7 

-29.3 


TABLE V. Reaction rates for the transition p p —> t ^Be 



calculated with the three different aN potentials |27-29| 


Potential [p^ 

m 

P [sec 

-1] 1688 1129 

238 
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Figure Captions. 

FIG.l. The correlation between the abundance of li and Fe in old stars (normalized to the 
abundance in the sun). 

FIG.2. The correlation between the abundance of Mg and Fe in old stars (normalized to 
the abundance in the sun). 

FIG.3. Yp denotes the predicted primordial abundance of ^He relative to hf as a function 
of T] for three different numbers of neutrino species Ni^. The other curves show the relative 
abundances of D, ^He, and Ti as a function of r]. 

FIG.4. New experimental! data (at energies less than 100/ceG) and old data (above 100/ceG) 
for the S-factor. The full calculation (solid line) is shown together with a calculation using 
only the {L = 2,S = 2) conhguration (short-dashed line). 

FIG.5. Jacobi coordinates for the system hZ/i. 

FIG.6. The spectrum of the ^Be nucleus (schematically) ^1 . 

FIG.7. Two different sets of Jacobi coordinates appropriate to describe the initial and the 
hnal states of the four-body system. 

FIG.8. a) The neutron is bound by the proton b) it is bound by the ^Be nucleus. 

FIG.9. The spectrum of the ^^Be nucleus (schematically) 0. 

FIG. 10. Transition through the transfer of two neutrons. 

FIG.11. Transition through the resonance state of the ^^B* nucleus. 

FIG. 12. The basic reactions of the muon catalyzed dt fusion cycle. 

FIG. 13. The excited state of ^Be and the position of the threshold for the decay into d+ 
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FIG. 14. The excited state of ^Be and the position of the threshold for the decay into d+ 
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FIG. 15. The basic reactions of the muon catalyzed d fusion cycle. 

FIG. 16. The excited state of and the position of the threshold for the decay into ^^-|-2p 


151 
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